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Abstract 

The weak radiative decay H° — > Ae + e~ has been detected for the first time. We 
find 412 candidates in the signal region, with an estimated background of 15 ± 5 
events. We determine the branching fraction B(E° — > Ae + e~) = [7.6 ± 0.4(stat) ± 
0.4(syst) ± 0.2(norm)] x 10~ 6 , consistent with an internal bremsstrahlung process, 
and the decay asymmetry parameter a=Aee = —0.8 ± 0.2, consistent with that of 
H° — > A7. The charge conjugate reaction H° — > Ae + e~ has also been observed. 
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1 Introduction 



Since the discovery of hyperons, their (weak) radiative decays have held par- 
ticular interest [1,2]. Still, the precise nature of the decays themselves remains 
an open question [3,4]. 
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Reliable techniques to predict branching ratios remain elusive. Furthermore, 
because SU(3) symmetry is broken only weakly in this regime, weak radiative 
decays should approximately conserve parity [5]. Consequently, the asymme- 
tries of decay angular distributions should be small. However, results from 
experiments indicate a relatively large (negative) asymmetry in every mode 
investigated [6]. A number of models have been proposed to explain this appar- 
ent discrepancy [7]. Experimental results tend to favor pole models or models 
based on chiral perturbation theory, which correctly find the sign of the asym- 
metry. Recently, a resolution of at least part of the puzzle has been offered [8] . 



When the NA48 Collaboration undertook investigations with a high-intensity 
Kg beam in 2002, trigger strategies for identifying radiative hyperon decays 
were included from the outset. The production over the full course of the run 
of more than 3 x 10 9 neutral cascades, S°(1315), offered NA48 unmatched 
sensitivity for the study of such decays £j 



This Letter details the measurement with these data of the weak radiative hy- 
peron decay 2° — > Ae + e~. This is the first measurement of this decay channel. 
If one assumes an inner bremsstrahlung-like mechanism producing the e + e~ 
pairs, the expected rate for this process may be estimated naively assuming 
the (virtual) photon converts internally (Dalitz decay) or by using the ma- 
chinery of QED as carried out in rate predictions for S° — ► Ae + e~ [1,2]. The 
results give a range from about 1/182 to 1/160 of the rate of S° — > A7, or 
(6.4 — 7.3) x 10~ 6 . Such a process should exhibit a decay asymmetry like that 
in 2° -> A7. 



2 Data 



The signal was sought among events containing one A (which decayed in-flight 
to a high-momentum proton and a much lower momentum 7T~), one electron, 
and one positron, all in time. An additional in-time photon was required for 
the normalization channel, (5° — ► Air , tt° — ► e + e~7), a relatively abundant 
process whose final state is similar to that of the signal and which was selected 
via the same trigger tree as the signal channel. 



1 The S° production rate was about 1/11 that of H . This Letter presents numerical 
results for the H° only. 
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2.1 Beam line and detector 



The NA48 beam line was designed to produce and transport both K\ and Kg 
beams simultaneously [9]. For the 2002 run, in order to increase dramatically 
the intensity of the K s beam, the K\ target was removed and the K\ beamline 
blocked, the proton flux on the K s target was greatly increased, and a 24 mm 
platinum absorber was placed after the Be target to reduce the photon flux 
in the neutral beam. An additional sweeping magnet was installed across the 
5.2-meter long collimator, which, tilted at 4.2 mrad relative to the incoming 
proton beam, selected a beam of long-lived neutral particles (7, n, K°, A, and 
S°). In each 4.8 s spill, occurring every 16.2 s, ~ 5 x 10 10 protons impinged 
on the target. Approximately 2 x 10 4 S°s, with momenta between 60 and 220 
GeV/c, decayed in the fiducial volume downstream of the collimator each spill. 

The detector for the 2002 run was identical to that used for NA48's measure- 
ment of direct CP- violation [9], except that the tagging counter immediately 
after the last collimator was removed. 

The neutral beam exited the final collimation into an evacuated tank, approx- 
imately 90 m in length, terminated by a Kevlar window 0.3% of a radiation 
length thick. The detector was arrayed immediately downstream of this win- 
dow. 

A magnetic spectrometer followed the decay volume. It consisted of four drift 
chambers, two before and two after an analyzing magnet which provided a 
transverse momentum kick of 265 MeV/c in the horizontal plane. The cham- 
bers were identical, with two planes of sense wires in each of four orienta- 
tions (x, y, u, v), vertical, horizontal, and at ±45°. The u and v wires of 
the third chamber were the only ones left un-instrumented. Track-time reso- 
lution was about 1.4 ns. Space-point resolution was approximately 150 /im 
in each projection, and the momentum resolution (with p in GeV/c) was 
&p/p = 0.48% © 0.015% x p. The resulting m n + w - resolution in K s — > 7r + 7r~ 
decays was 3 MeV/c 2 . 

A liquid krypton calorimeter (LKr) detected and measured the energy and 
position of electromagnetic showers. Its active region was divided transversely 
into approximately 2 cm x 2 cm cells, and its depth was 27 radiation lengths. 
Its single-shower time resolution was less than 300 ps; its transverse position 
resolution was better than 1.3 mm for a single photon of energy greater than 20 
GeV; and its energy resolution [10] was a(E)/E = 3.2%/ v / £©9%/£©0.42%, 
where E is in GeV. The resulting m 77 resolution in n° — > 77 decays was 
approximately 1 MeV/c 2 . 

The sensitive region of the electromagnetic calorimeter primarily constrained 
the fiducial volume of the experiment. Seven rings of scintillation counters 
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bounded, in projection, the edges of this acceptance region, and the last two 
rings acted as trigger vetoes of extraneous activity. 

A scintillator hodoscope, comprised of segmented horizontal and vertical strips 
arranged in four quadrants and located between the downstream end of the 
spectrometer and the upstream face of the calorimeter served as a zeroth-level 
charged-track trigger. Beyond the electromagnetic calorimeter stood an iron- 
scintillator sandwich hadron calorimeter and three layers of muon counters, 
each shielded by an iron wall. 

The entire detector array was sampled every 25 ns. An event trigger initiated 
a readout of information within a 200 ns window around the trigger time. In 
this way, time sidebands allowed investigations of accidental activity. 

The experiment employed a multi-level trigger designed to maximize flexibility 
while minimizing pile-up, dead-time losses, and the collection of uninteresting 
events. To be included in the present analysis, events passed the lowest level 
hardware trigger if a horizontal-vertical coincidence occurred in at least one 
quadrant of the scintillator hodoscope, there were no in-time hits in the veto 
rings, at least three views in the first drift chamber registered more than two 
hits (as required in the case of more than one track), and either the energy 
in the electromagnetic calorimeter exceeded 15 GeV or the total energy in 
the electromagnetic and hadron calorimeters exceeded 30 GeV. The next level 
trigger required more than one track to have passed through the spectrome- 
ter forming one or more good vertices^] The highest level trigger, an offline 
software cull, passed events in which a two-track invariant mass was consis- 
tent with that of a A and contained at least one high-energy cluster in the 
calorimeter not associated with either track forming the A. 

A downscaled sample of minimum bias events was collected concurrently with 
the physics data. Complete trigger information was available for these events, 
so trigger efficiencies could be measured. The relative fraction of events con- 
taining all signal final state particles that passed the required triggers was 
e ti?g = (96.5 ± 0.2)%, while the relative fraction of normalization events was 
e t n r ° g m = (97.1 ±0.2)%. 

2.2 Event selection criteria 

From events passing all trigger levels, those containing exactly four charged 
tracks, two of each charge sign, that passed well within the fiducial volumes 
of the first and fourth drift chambers were kept for further analysis. 

2 A good vertex is defined, in this context, as the occurrence of two tracks passing 
within 5 cm of one another between the target and the first drift chamber. 
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Signal event simulation showed that 99% of final state pions, electrons, and 
positrons had momenta of less than 30 GeV/ c. A track with momentum greater 
than 3 GeV/c and associated shower energy within 5% of this momentum was 
identified as an electron or positron, depending on charge. A positive track 
whose momentum was greater than 30 GeV/c and either had no associated 
electromagnetic shower or the shower energy to momentum ratio was less than 
0.8 was identified as a proton. If no such track was found, or if there were not 
both an electron and a positron identified, the event was abandoned. If the 
final track had a momentum greater than 4 GeV/c, but not more than 1/3.7 
that of the proton track, it was identified as a pion. Otherwise, the event was 
abandoned. 

The tracks associated with the proton and pion had to be separated by at least 
5 cm in the first drift chamber and their detection times had to be within 2 ns. 
If not, the event was abandoned. The distance-of-closest-approach (doca) of 
the two tracks when projected back towards the target was required to be less 
than 2.2 cm, and the longitudinal position of this doca had to lie between 4 
and 40 meters down stream of the target for the event to be further considered. 
The momentum vectors of the two tracks were projected, with respect to a 
reference frame centered on the beam axis, from their positions in the first 
drift chamber onto the face of the LKr. These projections were weighted by 
the relativistic energies of the particles associated with the respective tracks, 
added vectorially, and then normalized to the energy sum of the two particles. 
The result, a quantity called the center of gravity (COG), had to be greater 
than 8 cm to ensure that a parent of the two tracks was unlikely to have been 
directly produced in the target. The COG of a directly produced particle 
should be small. 

The invariant mass of surviving proton and pion candidate pairs was calcu- 
lated. If the result differed from the nominal mass of the A by more than 3 
MeV/c 2 (approximately 3a), the event was abandoned. 

The electron and positron tracks had to have times within 2 ns and a spatial 
separation in the first drift chamber of at least 2.5 cm. The latter require- 
ment rejects conversions in the Kevlar window. Any unassociated shower in 
the calorimeter with energy above 1.5 GeV disqualified the event as a signal 
candidate. 

A shower of between 3 and 120 GeV in the electromagnetic calorimeter was 
considered a photon candidate for the normalization channel if it was unasso- 
ciated with any track, centered within the fiducial volume of the detector at 
least 5 cm from a dead cell, and isolated from any other shower. 

Finally, for both signal and normalization channels, the event COG, which 
ideally would be (see above), had to be equal to or less than 6 cm. 
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A signal (normalization) region was defined as 2a either side of the nominal 
Ae + e~(7) invariant mass, where a m = 1 MeV/c 2 . For the A, the pir invariant 
mass was used. Selection from the entire data set according to these criteria re- 
sulted in 412 signal candidates and 29522 normalization events reconstructed. 



3 Acceptance and reconstruction efficiency 

The product of geometrical acceptance (A) and selection criteria efficiency (e) 
was determined with a Monte Carlo simulation. Nearly 10 5 signal-like events 
were generated according to a two-body model of a A and a virtual photon. 
The model included the decay parameter a = —0.78, found for the decay 
S° — > A7 [11], and a energy distribution for the converting photon, as 

would be the case for inner conversion. In this way, the model was intended 
to represent inner bremsstrahlung production. Generated events were stepped 
through a GEANT simulation of the NA48 detector and analyzed as real data, 
with the result: (A x e) s ; g = (2.69 ± 0.05)%. For the normalization channel, 
about 160 x 10 6 events (about 7x the measured flux) were generated with 
the latest PDG values for the decay parameters incorporated [6]. The result 
of the detector simulation and reconstruction was (A x e) norm = (0.1251 ± 
0.0003)%. Radiative corrections, using PHOTOS [12], were included, as was 
a H° polarization of —10% for signal generational 



4 Background 

Two sources of background were identified: physics and accidentally in-time 
combinations. 

4-1 Physics backgrounds 
4.1.1 S ^Att 

The S° decays predominantly to Air . If the 7r° Dalitz-decays, and the photon 
goes undetected, the final state is that of the signal. Similarly, if the 7r° decays 
via the double-Dalitz mechanism, and an electron and a positron go unde- 
tected, the final state is again that of the signal. Finally, the tt° — > e + e~ decay 
results in an irreducible background, but its rate is very small. Simulations 

3 This polarization value is consistent with that reported by other experiments [13] 
and with indications from an ongoing study of the NA48 beam. 
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of each of these channels at about seven times the flux lead to estimates of 
4.6 ± 0.8, 0.1 ± 0.1, and 1.2 ± 0.4 events, respectively, infiltrating the signal 
region. 



4-1.2 Kaon decays 

The flux of neutral kaons was an order of magnitude larger than that of the 
S°. The decay Kg — > 7r + 7r~e + e~ has a branching fraction of 4.7 x 10~ 5 . If 
one of the pions met the requirements of a proton in this analysis, and the 
resulting; then this process would mimic the signal. Simulation with 

twice the flux of such events demonstrated that an explicit mass cut | m nnee — 
m K o |> 0.015 GeV/c 2 eliminated essentially any trace of this background with 
negligible impact on signal-finding efficiency The decay chain K\ — > 7r + 7r~7r°, 
7T° — > e + e~7, has a product branching ratio of about 1.5 x 10~ 3 . The K® L 
lifetime and their typical momentum of 80 GeV/c mean that about 4% of them 
decay in the experiment's decay volume. For these to become a background to 
the Aee signal, a pion would have to be mistaken as a proton and the invariant 
mass of it combined with that of the other pion would have to be close to that 
of the A. In addition, the photon would have to go undetected. Because of 
this last condition, an explicit kaon mass cut would be ineffective in reducing 
the background. On the other hand, the efficiency for this chain appearing in 
the signal region is correspondingly reduced and the COG is smeared out. We 
estimate on the basis of Monte Carlo simulation that 2 ± 2 such events will 
populate the signal region. 



4-1.3 Accidentally in-time combinations 

We estimated the contamination by accidental coincidences four ways: 

(1) Running the same analysis on the data, but requiring that the final-state 
leptons have the same charge. 

(2) Requiring that at least one track or shower be between 10 and 20 ns 
out-of-time and scaling appropriately. 

(3) Taking events with m pn values between 7 and 10 standard deviations 
from the central value (ttia) and computing m\ ee . 

(4) Defining two "side-band" regions, one along each axis in COG-versus- 
mA ee space [see, in Figure 1, the hatched rectangles at high COG and 
high mass; each region has the same "area" as the signal region, the open 
rectangle in the figure]. 

These approaches, which are not independent, yielded between 1 and 9 events 
in the signal region; we take the number to be 7 ± 5 events. 
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Fig. 1. COG versus m\ ee after all other selection criteria were imposed. The three 
hatched boxes are side-band regions. The signal region is the open box at low COG 
around m=o. The side-band regions at high mass-low COG and high COG were used 
to estimate accidental and combinatoric backgrounds in the signal region. All three 
side-band regions were used in the subtraction of background under the decay-angle 
distribution (see text). 



In conclusion, combining the physics backgrounds with those attributed to 
accidentals and combinatorics, the estimated number of background events in 
the signal region is 15 ± 5 [see Table 1 for a summary of the background 
estimation] . 



The background contamination of the normalization sample was estimated 
from the tails of the m eei spectrum, which peaks sharply at m n o. Including 
a linear extrapolation under the mass peak, the number was estimated to be 
428 ± 258. 
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Table 1 



Sources of expected background events. 



Source 


Estimate 


~° Avr°, tt° e+e- 7 


4.6 ±0.8 


~° Avr°, ^° e+e-e+e- 


0.1 ±0.1 


~o ^ A7r0) ^ e+e - 


1.2 ±0.4 


Kaon Decays 


2 ±2 


Accidentals & Combinatorics 


7±5 


TOTAL 


15 ± 5 



5 2° Flux 



The total number of S° produced during the run was estimated by fully re- 
constructing S° — > A7T°, 7T° — > e + e _ 7 events without a longitudinal vertex 
position cut and using the equation 

-r N norm N norm bkgd / 1 \ 

H ° = (A x e) norm e^ m £(S° -+ Att0)S(A -+ p7r-)£(7r° -+ e+ e - 7 ) 1 j 

From the entire data set, 29522 such events were reconstructed. After back- 
ground subtraction, this gives an integrated flux of 

$ H o = (3.15 ± 0.03 ± 0.08) x 10 9 . 



The first uncertainty is due to statistics, and the second is from branching 
fraction uncertainties, primarily that on B(n° — > e + e~ r y). 

Table 2 

Quantities that entered into H° flux calculations. 



No. of events in signal region 


29552 


Estimated no. of background events 


428 ± 258 


[A ~X. e)norm 


(0.1251 ±0.0003)% 


norm 
trig 


(97.1 ±0.2)% 


B(E° Avr°) 


0.9952 ± 0.0003 




0.639 ±0.005 


B(tt° e + e- 7 ) 


0.01198 ±0.00032 
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6 Results 



At the end of the analysis, 412 events were found in the signal region [see 
Figure 2]. 




1.29 1.3 1.31 1.32 1.33 1.34 



m Aee (GeV/c 2 ) 

Fig. 2. mAee after all selection criteria. Arrows indicate signal region. Stacked in 
various hatchings (see legend) are the estimated sources of background. 
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Fig. 3. Reconstructed m ee spectra from data (points), (solid line), 2-body 

flat (dashed line), and 3-body phase space (dotted line). The distributions from 
simulated data have been normalized to contain the same number of events in the 
signal region of the data, without background subtraction. 



6.1 m ee spectrum 



The associated m ee distribution is consistent with a 1/m^ shape [see Fig- 
ure 3] , and we consider only this model (presumably inner bremsstrahlung) in 
determining of the branching fraction, including systematic uncertainties. 



14 



6.2 Branching fraction 



Given the background estimate, efficiencies, and flux discussed above, and 
the PDG A — > pir~ branching ratio [see Table 3], the branching ratio for 
2° — > Ae + e~ is determined to be 

B(E° -> Ae+e-) = (7.6 ± 0.4) x 1(T 6 

where the uncertainty here is statistical only. 
Table 3 



Quantities that entered into branching fraction calculations. 



No. of events in signal region 


412 


Estimated no. of background events 


15 ±5 


{A x e) sig 


(2.69 ±0.05)% 


sig 
trig 


(96.5 ±0.2)% 


2° flux 


(3.14 ±0.03) x 10 9 


£>(A — > pir~) 


0.639 ±0.005 



A check for S° — > Ae + e~ found a clear peak of 24 events, of which, roughly, 
as many as 7 may be background, in the invariant mass plot. This number is 
consistent with the S° branching fraction and relative S° and S° production 
rates. The production mechanism, kinematics and backgrounds of the H°, 
however, differ from those of the S°, and no further consideration of this 
charge-conjugate channel is given here. 

6.3 Systematic uncertainties 

Analysis selection criteria were varied when looking at the data and when 
determining reconstruction efficiencies. The branching fraction result was most 
sensitive to the treatment of the reconstructed S° vertex and backgrounds from 
the S° — > A7T° channel in relation to m ee . No cut was placed on the longitudinal 
position of the S° vertex. Requirements varying the minimum longitudinal 
position of the vertex in 6-m intervals beginning before the target (to account 
for resolution effects) resulted in branching fraction changes of between 0.2% 
and 3%. We assign the highest variation (±3%) as a systematic error. 

It was possible to eliminate nearly all physics backgrounds by excluding signal 
events with 0.100 GeV/c 2 < m ee < 0.135 GeV/c 2 , which, according to signal 
Monte Carlo, reduces the reconstruction efficiency by 5%. Cutting this region 
from the final data sample, and recalculating the branching ratio, results in a 
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shift of 1.8%, which was included symmetrically as a systematic uncertainty. 
These, along with smaller variations in the branching fraction resulting from 
other modifications of the selection criteria, were added in quadrature to give 
a systematic uncertainty of ±3.6% on the branching fraction. 

We conservatively assign a relative ±1% uncertainty on the determination of 
the background to account for correlations in methods for estimating acciden- 
tally in-time events. 

The branching fraction differed by about 1% when signal and normalization 
modes were simulated with and without radiative corrections, and we include 
this difference symmetrically as a systematic uncertainty. 

For the Axe determinations, the S° polarization of simulated events was set 
to —10%. Samples of simulated data, generated with the polarization varied 
between 0% and —20% (±10%), were used to recalculate the branching frac- 
tion vary. The largest variation among these trials was 2.7%, and this variation 
is taken symmetrically as a systematic uncertainty. 

The decay asymmetry used in generating simulated signal events was that of 
the process S° — > A7 [11]. Our measurement, discussed below, is in agreement 
with this value, but with a 25% uncertainty. Varying our simulation within 
this 25% range changed the branching fraction by at most 2.5%, and this is 
symmetrically assigned to systematic uncertainty. 

The determination of the trigger efficiency and S° flux were discussed above. 
The difference between trigger efficiencies for signal and normalization chan- 
nels is taken as an uncertainty, affecting the branching ratio by 0.6%. An 
alternative, less direct, calculation of the flux was statistically consistent with 
the one described above. The two differed by 1.9%, and we conservatively 
include, symmetrically, this amount as a systematic uncertainty. 

The total systematic uncertainty on the branching fraction, recounted in Ta- 
ble 4, is ±5.7%, the sum in quadrature of each of the sources described. This 
gives a final branching fraction of: 



B(E° -> Ae + e-) = [7.6 ± 0.4(stat) ± 0.4(syst) ± 0.2(norm)] x 10~ 6 . 
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Table 4 

Sources of systematic uncertainty on the branching fraction. 



Source 


Fractional 




Uncertainty 


Detector Acceptance 


3.6% 


Background 


1.0% 


Radiative Corrections 


1.0% 


Polarization 


2.7% 


Signal Modeling 


2.5% 


Trigger Efficiency 


0.6% 


~° Flux 


1.9% 


TOTAL 


5.7% 



6.4 Asymmetry parameter 



The angular distribution of the proton relative to the S° line of flight in the 
A rest frame is given by [6]: 

dN N 

^^ = y (l-aHAeea_cosM- (2) 



The cos Op- spectrum from signal events was corrected by subtracting scaled 
backgrounds from the side-band regions indicated in Figure 1 and by dividing, 
bin-by-bin, the acceptance as determined from a 2° — > Ae + e~ simulation 
where the spectrum was generated to be flat in cos 9 p w. A two-parameter fit to 
this corrected spectrum gives the product of asymmetry parameters a^\ ee a_, 
where a_ is the asymmetry parameter for the decay A — > piT~ . This latter was 
taken to be a_ = 0.642±0.013 [6]. The fit (over the interval -0.8 < cos 9 pS < 1) 
[see Figure 4] to the data yields, 

a S Aee = -0.8 ± 0.2 



This is consistent with the latest published value of a=A 7 = — 0.78±0.18(stat)± 
0.06(syst) [11]. 
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Fig. 4. Background-subtracted and acceptance-corrected cos 9 P ~, distribution. The 
line is the fit result. 



7 Summary and conclusions 



The weak radiative decay channel S° — > Ae + e has been identified. Its branch- 
ing fraction has been determined to be 

B(E° -> Ae+e-) = [7.6 ± 0.4(stat) ± 0.4(syst) ± 0.2(norm)] x KT 6 , 
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consistent with an inner bremsstrahlung-like production mechanism for the 
e + e~ pair. The consistency is further supported by the m ee spectrum. The 
decay parameter 

«HAee = "0.8 ± 0.2, 

is consistent with that measured for S° — > A7. 

The existence of the charge conjugate reaction S° — > Ae + e~, has been con- 
firmed. 
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